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ABSTRACT 

We combine VI photometry from OGLE-III with VVV and 2MASS mea- 
surements of E{J — Kg) to resolve the longstanding problem of the non-standard 
optical extinction toward the Galactic bulge. We show that the extinction is 
well-fit by the relation Ai = 0.7AQ5xE{V - /) + 1.3700xE(J - Ks), or, equiva- 
lently, Aj = 1.217xE{V - + lA2Qx{E{J - Ks)/ E{V - I) - 0M33)). The 
optical and near-IR reddening law toward the inner Galaxy approximately fol- 
lows an Rv ~ 2.5 extinction curve with a dispersion a^y ~ 0.2, consistent with 
extragalactic investigations of the hosts of type la SNe. Differential reddening is 
shown to be significant on scales as small as as our mean field size of 6', with the 
1(7 dispersion in reddening averaging 9% of total reddening for our fields. The in- 
trinsic luminosity parameters of the Galactic bulge red clump (RC) are derived to 
be {Mi^Rc,(Ti,RC,o,{V-I)RC,o,<Jiv-i)ncdJ-Ks)RCfi) = (-0.12,0.09,1.06,0.121, 
0.66). Our measurements of the RC brightness, brightness dispersion and num- 
ber counts allow us to estimate several Galactic bulge structural parameters. We 
estimate a distance to the Galactic center of 8.20 kpc, resolving previous discrep- 
ancies in distance determinations to the bulge based on 7-band observations. We 
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measure an upper bound on the tilt a ~ 40° between the bar's major axis and 
the Sun-Galactic center line of sight, though our brightness peaks are consistent 
with predictions of an N-body model oriented at a ~ 25°. The number of RC 
stars suggests a total stellar mass for the Galactic bulge of 2.OxlO^°M0, if one 
assumes a Salpeter IMF. 

Subject headings: Galaxy: Bulge, fundamental parameters, stellar content, struc- 
ture - ISM: dust, extinction 



Introduction 



The central bulge of the Milky Way Galaxy is the only stellar spheroid for which we can 
measure detailed abundances, ages and all six phase space dimensions for individual stars, 
as well as the luminosity function and spatial distribution for the population as a whole. 
Some ~10% of the Milky Way's stars are bulge stars, including a disproportionate number 
of the oldest and most metal-rich stars. It is therefore evident that any theory of Galaxy 
formation and evolution is required to reproduce the observed properties of the bulge, and 
conversely, that the properties of the bulge should be measured as precisely and accurately 
as possible to best discriminate between different Galaxy formation models. 

However, as scientifically desirable as this greater project may be, it is also difficult 
due to several challenges that prevent further, deeper understanding of the bulge. There are 



significant correlations between chemistry and kinematics (Babusiaux et al. 2010 Uttenthaler 



et al.||20l"2 Ness et aL]|2012 ), distinct chemical subgroups ( Andrews et al.||2012 ), gradients in 



metallicity (Zoccali et al. 2008 ), deviations from the classical picture of the triaxial ellipsoid at 



large separations from the minor axis ( Benjamin et al.||2005 Cabrera-Lavers et al.|2008| ), and 
both large QNataf et al.]|2010| [McWiUiam fc Zoccah|[20lo| and small ( |Alard||2001t [Nishiyama 



et al. 2005) separations from the plane. These issues necessitate larger data sets and better 



models. The viewing angle a between the bar's major axis and the sun-Galactic center line of 
sight remains undetermined, with best-fit values ranging from from a = 13° ( Cabrera-Lavers 



et al. 2007) to a = 44° (Benjamin et al. 2005). This prevents further understanding of the 



inner Galaxy's gravitational potential, as uncertainties in the value of a are degenerate with 



those of the bar's axis ratios (Stanek et al. 1997) and rotation speed (Shen et al. 2010). 



^ Based on observations obtained with the 1.3 m Warsaw telescope at the Las Campanas Observatory of 
the Carnegie Institution for Science. 
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The most significant source of uncertainty, however, is the extinction. It averages 



3 toward the Galactic center (Schodel et ah 2010), suggesting Ay ~ 50 (Nishiyama et aL 



2008). For most of the bulge, values of Ay = 2 are typical (Sumi 2004). The high values 



of reddening close to the plane render it quite difficult to obtain spectroscopic observations, 
proper motions, and stellar density maps. Further from the plane, these can be obtained, 
but not fully understood due to significant zero-point uncertainties in the extinction, and 
indirectly, the distance. 

A further complication arises from the fact the extinction toward the inner Galaxy is not 



only large but also non-standard. This was first suggested by Popowski (2000) as a solution 



to the anomalous colors of bulge RR Lyrae (Stutz et al. 1999) and red clump (RC) stars 



(Paczynski et al. 1999 ). Gould et al. (2001 ) and Udalski (2003b ) were the first to demonstrate 



that the reddening law toward the inner Galaxy is described by smaller total-to-selective 
ratios than the "standard" values measured for the local interstellar medium, implying a 



steeper extinction curve and thus a smaller characteristic size for dust grains (Draine 2003). 



Udalski] ( |2003b[ ) measured dAj/dE{V - J) ^ 1.1 (denoted AAi/AE{V - I) in that work) 
toward several bulge fields, much smaller than the value of dAj/dE{V — /) ?a 1.45 suggested 
by the standard interstellar extinction curve of Ry = 3.1 (Cardelli et al. 1989 O'Donnell 



1994). Udalski (2003b) showed that applying the same methodology to observations of 



the Large Magellanic Cloud taken with the same instruments yielded dAj/dE{V — I) 
1.44, the standard value, demonstrating the robustness of the result. Further, not only 
was the reddening law toward the bulge found to be non-standard, it was also found to be 
rapidly varying between sightlines, with values of dAj/dE{V — I) ranging from 0.94 ± 0.01 
to 1.16 ± 0.03. The steeper extinction law toward the inner Galaxy has been subsequently 



confirmed with observations using Hubble Space Telescope (HST) optical filters by Revnivtsev 
et al. (2010), by analysis of RR Lyrae stars in OGLE-III ( Pietrukowicz et aL]|2012 ), and also 



in the near-IR (Nishiyama et al. 2008 Goshng et al. 2009 Schodel et al. 2010). Meanwhile 



Zasowski et al. (2009) and Fritz et al. (2011) both found that the extinction law toward the 



inner Galaxy was shallower (greyer) in the mid-IR. 

The variations in both the extinction and the extinction law made it difficult to reliably 



trace the spatial structure of the bulge (Majaess 2010). Applying the VI extinction maps 



of Sumi (2004) to the bulge color- magnitude diagram (CMD) implied a distance to the 
Galactic center of ~9 kpc (Rattenbury et al. 2007a VanhoUebeke et al. 2009), a large value 



relative to the geometrically determined distances to the Galactic center of 7.94 ± 0.42 kpc 



(Eisenhauer et al. 2003), 8.27 ± 0.29 kpc (Schonrich 2012), and 8.4 ± 0.4 kpc (Ghez et al. 



2008). As the structure of the inner Galaxy is a very sensitive probe of the environmental 



conditions in which the Galaxy formed and evolved (Athanassoula 2005 Inoue & Saitoh 



2012), an accurate spatial determination of the bar's morphology would yield a powerful 
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test of Galaxy formation models. Moreover, investigations of the metallicity distribution 
function of bulge giants has had to rely on imprecise and potentially inaccurate estimates of 
surface gravity and photometric temperature, furthering reducing our ability to probe the 
primordial conditions of the Galaxy. 

We resolve these issues in this investigation by combining OGLE-III observations in 



V and I with VVV and 2MASS measurements of E{J - K^) ( [Gonzalez et al.||2012[ ). We 
confirm previous findings that the VI extinction toward the inner Galaxy is steeper than 
standard, but also show that it is a little less steep than previously assumed. We show 
that this is likely due to an effect we label "composite extinction bias", which makes it 
unphysical to extrapolate a slope of dAj/dE{V — I) to E(y — /) = 0. Our parameterization 
for the extinction, Aj = 0.7465xE{V — I) + 1.3700xE{J — Kg), is less sensitive to the rapid 
variations in the extinction law than the computation of slopes dAi/dE{V — I). Whereas 
the latter must be computed from an ensemble of measurements spread across 30' or more, 
the former can be directly measured for each ~ 6' x 6' sightline. 

The structure of this paper is as follows. We summarize the data used in Section [2j 
Our methodology for measuring the parameters of the RC is described in Section [3} and we 
derive the intrinsic RC luminosity parameters in Section |4} We briefly state the properties 
of reddening that would be expected using a standard reddening curve in Section [5] The 
reddening measurements are presented and discussed in Section [6| including comparisons 



to the reddening maps of Schlegel, Finkbeiner & Davis (SFD, Schlegel et al. 1998) and the 
derivation of an empirical rule to estimate differential reddening. Methods to convert the 
reddening into an extinction using a single color are demonstrated to inevitably fail in Section 
[7} and a more successful extinction law is derived in Section [8] by including information 
from both E{V — I) and E{J — Kg). We demonstrate that reddening constraints from 



MACHO photometry have been misinterpreted in Section [9j In section 10, we show that our 
dereddened apparent magnitudes suggest a distance to the Galactic center Rq = 8.20 kpc, 
and a tilt between the Galactic bar's major axis and the sun-Galactic center line of sight no 



greater than a ~ 41°. We analyze our number counts for the RC in Section 11 and show 
that combining these with the assumptions of standard stellar evolution and a Salpeter IMF 
yields an estimated Galactic bulge stellar mass of M = 2.OxlO^°M0. The thickness of the 



Galactic bulge is discussed in Section [T2j The data structure is summarized in Section [13 
Results are discussed in Section 



- 5 - 



Data 



OGLE-III observations were taken with the 1.3 meter Warsaw Telescope, located at the 
Las Campanas Observatory. The camera has eight 2048x4096 detectors, with a combined 
field of view of 0.6° x 0.6° yielding a scale of approximately 0.26"/pixel. We use observations 
from 263 of the 267 OGLE-III fields directed toward the Galactic Bulge, which are almost 
entirely within the range —10° < / < 10° and 2° < |6| < 7°. We do not use 4 of the fields, 
BLG200, 201, 202, and 203; located toward (/, 6) ^ (-11°, -3.5°), due to the much higher 
differential reddening and disk contamination toward those sightlines. The photometric 
coverage used in this work is shown in Figure [1} Of the 263 fields used, 37 are toward northern 
latitudes. More detailed descriptions of the instrumentation, photometric reductions and 



Udalski 


(2003a 


), 


Udalski et al. 


(2008 


) and 


Szymahski 



et al. (2011). OGLE-III photometry is available for download from the OGLE webpage]^ 




Fig. 1. — Coverage of the OGLE-III Galactic bulge photometric survey used in this work, 
overplotted on an optical image of the same area. Galactic coordinate system shown. Red 
squares denote OGLE-III fields used in this work, and yellow squares denote fields not used. 



We also make use of data from the Two Micron All Sky Survey (2MASS, Skrutskie 



^http://ogle. astrouw.edu. pi/ 
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et al. 2006), which we cahbrate on the measurements of Gonzalez et aL (2012), who used 
data from the VISTA Variables in The Via Lactea survey (VVV, [Saito"eraL][20l2| ). The 
calibration is discussed in Section l3^ 



3. Measuring the Red Clump 

The RC is a prominent, well-populated, and localized feature of Galactic bulge CMDs 
(Terndrup 1988 Stanek et al.||1994 ), of which we show two examples in Figure [2] The color. 



color- dispersion, apparent magnitude, magnitude dispersion and normalization of the RC 
vary from sightline to sightline, rendering it a sensitive probe of the reddening toward the 
bulge, its distance, and its underlying structure. We measure these properties across the 
OGLE-III bulge sky in the following manner. 

Each of 2,104 OGLE-III subfields (eight detectors over 263 fields) used in this work was 
split into 1, 2, 3, 6, 8, 10, 15 or 21 rectangles. Smaller rectangles were used toward regions of 
the sky where the surface density of stars was higher. Our average rectangle size is ~ 6' x 6'. 
The total number of sightlines is 9,744, though for most of this work we only make use of 
9,014 of those sightlines that are no nearer than 7' or 3 half-light radii to a known Galactic 
globular cluster]^ that satisfy our photometric completeness criteria of (^ — /)i?c < 3.30 and 
Irc < 17.70, and that are not flagged as being highly differentially reddened, or otherwise 
problematic. We publish the best-fit parameter values for the remaining sightlines but do 
not incorporate them in our analysis. 

The complexity of bulge CMDs requires that we be careful before fitting a luminosity 
function. The typical limits to the color-magnitude selection box are given by: 

- 0.30 < (F - /) - (F - /)«c (1) 

-1.5</-/^,c<1.5. (2) 

The purpose of the color-magnitude selection is to select as many bulge giants as possible 
without selecting too many foreground stars, which populate a sequence of stars ~0.6 mag 
bluer than the RC at the brightness of the RC, but that merges with the bulge red giant 



(RG) branch at a luminosity ~2 mag fainter than the RC (Kiraga et al. 1997). The color 
range is therefore adjusted when the foreground main-sequence stars either have very similar 
or very distinct colors from the bulge RG branch. 



■^As searched for using A Galactic Globular Cluster Database: 
which is based on the Harris catalog ( Harris|[l996, 2010 edition I. 



http: / /gclusters. altervista.org/index.php 
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3.1. Mean Magnitude, Magnitude Dispersion and Normalization 

of the Red Clump 

The standard methodology for investigations of the RC toward various stellar systems 



is the Paczynski-Stanek equation (jPaczynski & Stanek||1998|): 

N{I)dI = a + b{I 



+ c{I - / 



RC) 



+ 



RC 



exp 



2 



IRC) 



24c 



(3) 



where N{I)dI is the number of stars at magnitude I in an interval of length dl, the Gaussian 
parameters Irc, ctrc, and Nrc measure the mean magnitude, magnitude dispersion, and 
number of RC stars, and a quadratic polynomial is fit for the underlying luminosity function 
of red giant (RG) stars. Though the broad application of the Paczynski-Stanek equation 
demonstrates its versatility, we modify it to enhance our accuracy and precision. 



As in some previous works (Nataf et al. 2010, 2011a|[b|[c Nataf & Udalski 2011), we 
fit the luminosity function of the RG branch to a 2-parameter exponential rather than a 
3-parameter quadratic. The reduced number of free parameters makes the fitting routine 
more stable. We found that fitting a quadratic to the RG branch can lead to catastrophic 
errors. That is because for large values of aRc, the Gaussian becomes degenerate with the 
quadratic term, leading to even larger values of aRc and Nrc at the expense of an unphysical, 
negative normalization to the RG branch. The exponential satisfies the physically-motivated 
condition of being both a strictly positive and strictly increasing function of magnitude, which 
stabilizes it. It is also sound theoretically, as stellar models actually predict an exponential 
distribution to the magnitudes of RG stars outside the red giant branch bump (RGBB) 



(Castellani et al. 1989). We also accounted for the RGBB and asymptotic giant branch 



bump (AGBB). We parameterized the luminosity function as follows: 



N{I)dI = A exp 



IRC 



N 



RC 



N^ 



RGBB 



'^ttctrgbb 



exp 



RGBB 



RGBB 



+ 

Nagbb 



exp 



RC) 



24c 



+ 



2vr(TACBB 



exp 



{I — IagbbY 

AGBB 



(4) 



We fixed the parameters of the RGBB and AGBB to their mean values (from Nataf et al. 



2011c) to minimize the systematic effect of extra parameters on the fits: 

Nrgbb = 0.201 X Nrc 
Nagbb = 0.028 x Nrc 
I RGBB = Irc + 0.737 
I AGBB = Irc — 1-07 

O'rgBB = AGBB = o'rc 



(5) 
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We imposed the constraint that the integral J N{I)dI be equal to the number of stars in 
each fit, and thus we have four free parameters: Irc, o'rc: B, and Nrc/A. We contrast 
two CMDs in Figure |2] specifically selected to demonstrate the effects of reddening and 
extinction. In Figure |3} we show the CMD for a typical sightline, as well as our color- 
magnitude selection box for the fit, the corresponding magnitude histogram and best-fit to 
the RC+RG+RGBB+AGBB luminosity function in /. 
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Fig. 2. — OGLE-III CMDs for two distinct sightlines shown on the same figure. The morphol- 
ogy of the CMDs are similar, with both having a foreground disk main-sequence component 
to their left with a bulge RG branch including a RC to their right. However, the stars 
toward (/, b) = (2.14°, —1.89°) (red), are redder and fainter than the analogous stars toward 
(/,6) = (2.14°, —4.10°) (blue), due to higher interstellar extinction. The RC toward the 
former sightline is measured to be 1.17 mag redder and 1.44 mag fainter. 
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1.5 2 2.5 3 



OGLE V-I 



Fig. 3.— OGLE-III CMD toward {l,b) = (-2.29°, -3.12°) shown in the left panel. The 
best-fit values of the color and magnitude of the RC, the magnitude dispersion and the 
color dispersion are shown on the top-left of the left panel. The color-magnitude selection is 
denoted by the thick black lines. The magnitude histogram of stars in the color-magnitude 
selection box is shown to the right of the CMD, on the same scale to the vertical axis, along 
with a model fit for the RG, RC+RG and total RG+RGBB+RC+AGBB. 
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3.2. {V — I) Color and Color-Dispersion of the Red Clump 



The color of the RC is first assumed to be that of the nearest of 92,000 values of iV—I)Rc 



measured for the study of Nataf et al. (2010 ). The measurement of {V — I)rc is then slightly 



adjusted in this work. We start with the old value, and then measure the position of the RC 
in / by fitting Equation |4j Two colors are subsequently fit for, one for the RC, and one for 
the position of the foreground disk component at the brightness of the RC (typically ~0.6 
mag bluer than the RC), easily discernible in Figures |2] and [s] Every star in the CMD is 
then assigned to the closer of these two colors, and the two colors selected are those that 
minimize the weighted variance of the difference between the colors of stars and the nearest 
of the two trial colors (RC and foreground disk component), where the weights Wi are given 
by: 



Nrc 



exp 



A exp 



B{I - 1 



RC) 



Nrc 



exp 



(6) 



and thus in practice only stars satisfying |/ — I^d ^ o"/?c contribute to the fit to the color 
and color-dispersion. We recursively removed 2.5a outliers from the fit to {V — I)rc- The 
values of the color {V — I)rc and color-dispersion a(^v-i)Rc adopted once no outliers 
remain. If the resulting color disagrees with the color first assumed by 0.03 mag or more, 
we redo the fit to the magnitude distribution and subsequently to the color with an updated 
color-magnitude selection box. If the fit fails to converge after 5 iterations, the 0.03 mag 
condition is relaxed, and the sightline is flagged as problematic and not incorporated within 
any of our analysis. 



3.3. (J - K,)rc of the Red Clump 



Gonzalez et al.| ( |2012 ) measured (J — Ks)jic across the viewing area (—10 < / < 



-1-10.2, —10 < b < +5) using data from the VVV survey (Saito et al. 2012). These results 



are used in our work, and cover most of the OGLE-III bulge viewing area. 

For the remaining part of the sky, we used {J — Ks)rc measurements from 2MASS data 



(Skrutskie et al. 2006) that we calibrated on the measurements of Gonzalez et al. (2012). 



The magnitude limit of 2MASS is typically brighter than the bulge RC. However, 2MASS 
does reach the RG stars that are at the color of the RC and brighter than the RC. We cross- 
matched the OGLE-III and 2MASS source catalogs, yielding a VIJKg photometry list for 
all RG stars brighter than the RC. For RG stars satisfying \ (y — I)rg — {V — I)rc\ < 0.33, 
we computed the linear relations {I — J),{I — Kg) = a + bx ({V — I)rg — (V — I)rc)- The 



difference of the two relations then yields (J — Kg) for RG stars at the {V — I) color of the 
RC. This methodology was first developed as a means to study gravitational microlensing 
events toward the bulge (Bennett et al. 2010; Gould et al. 2009). 



0.15 I r 
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(v-i) 



RC.OGLE-III 



Fig. 4. — Difference in {J — Ks)rc inferred from 2MASS photometry with that directly 
measured from deep VVV photometry by Gonzalez et al. (2012). The relationship is (J — 
Ks)kc,2MAss -{J - Ks)^cyvw = -0.035 + 0.031 x (l^ - /)rc,ogle-iii, which we subsequently 
correct for. Points denote each measurement. Both the medians to 50 bins and the best-fit 
line are shown by thick black curves. 



4. Calibration of the Intrinsic Luminosity Parameters of the Red Clump, 

{V - I)rc\o, Mi^rc, (^i,rc,o and (J - Ks)RCfl 

We adopt {Mi^Rc,ai,Rc,o.{V-I)Rc,o.(^{v-i)nc^iJ-Ks)RCfi) = (-0.12,0.09,1.06,0.121, 
0.66) for the mean absolute magnitude and magnitude dispersion in J, intrinsic {V — I) color, 
intrinsic (V" — /) color-dispersion, and intrinsic (J — Kg) color of the RC. Our justifications 
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for Mj^Fic, {V — I)RCfi, and (J — Ks)RCfl is provided below. The intrinsic color dispersion 
of the RC is measured to be a(v-i)iico ~ 0-121 in Section [7| and an estimate of the intrinsic 
magnitude dispersion ai^RCfi ~ 0.14 — 0.17 was provided by Nataf & Udalski (2011), but is 
revised downward to o"/,kc,o = 0.09 mag in this work. 



4.1. {V — I)rc,o From the Spectroscopy of Microlensed Dwarf and Subgiant 

Stars 



Spectroscopic investigations of high-magnification microlensing events toward the bulge 
have yielded detailed abundances and reddening-independent temperatures for bulge main- 



sequence turnoff and subgiant stars (Johnson et al. 2007, 2008 Cohen et al. 2008, 2009 



Epstein et aL]|2010t [Bensby et ar]|2009| |2010| |201lD . The intrinsic {V - I) colors of these 
stars can be estimated from the measured temperatures and metallicities and the empirical 



[Fe/H]-color-Teff calibration of Casagrande et al. (2010 ). This intrinsic color can be compared 
to the observed color to obtain the reddening. The assumption that the reddening to the 
source equals the reddening to the RC centroid yields {V — I)rc,o = 1-06 (Bensby et al. 
20TT|. 



4.2. {V - I)rc\o and Mj^rc by Linear Interpolation with 47 Tuc and NGC 6791 



We measure the properties of the RC for two old stellar populations that are well- 
calibrated to provide an empirical relation. 

We first use photometry for the Galactic globular cluster 47 Tuc (NGC 104) from 

the catalog of Sarajedini et al. (2007). The VI magnitudes are obtained by convert- 



ing from the original F606W and F81AW photometry taken with the (HST). The limits 
used to identify RC stars are the same as in Nataf et al. (2011b): (0.82 < V — I < 
0.98, 13.80 < < 14.20), thus delineating a box on the CMD containing 546 stars. We 



measure {{V -I)rc, Irc, ct/,rc)47Tuc = (0.913, 13.090, 0.062). Thompson et al. (2010) pre 



cisely measured the properties of an eclipsing binary (EB) pair in the cluster, and determined 
(m - M)v = 13.35 ± 0.08. Assuming a cluster reddening E{B - V) = 0.04 (|Harris|[T996y 



2010 edition) and a standard reddening law (Cardelli et al. 1989 [O'Donnell 1994), we find 
{{V - I)RCfi, Mi^Rc, a7,Rc,o)47 Tuc = (0.861, -0.208, 0.062). For the cluster metalhcity, we 



adopt [Fe/H]= -0.76 (Koch fc McWilham 2008). 



For the open cluster NGC 6791, we use the photometric catalog of Stetson et al. (2003), 



which was updated for the work of Brogaard et al. (2012) and now includes corrections for 
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differential reddening. We draw a box around the RC using the hmits (1.25 < V — I < 
1.40, 13.05 < / < 13.35), thus dehneating a box on the CMD containing 22 stars. We 



measure {{V - I)rc, Irc, o-/,i?c)NGC 679i = (1-326, 13.243, 0.050). Brogaard et ah (2011 



2012) precisely measured the properties of 3 eclipsing binary pairs in the cluster as well as 
the cluster CMD, and found EiV-I) = 0.174, {m-M)v = 13.51±0.06, and [Fe/H]= +0.29. 
We thus iinfer ((V^ - J)^c,o, M,,^c, c^/,iJC,o)NGC 679i = (1-152, -0.093, 0.050). 

The assumption of linear population effects in the metallicity range bracketed by 47 Tuc 
and NGC 6791 yields the following empirical relations: 

{V - I)b,c,o = 1-09 + 0.277 x ([Fe/H] - 0.05) (7) 
Mj^RC = -0.12 + 0.110 X ([Fe/H] - 0.05) (8) 



The bulge RC has a measured metallicity distribution function with mean [Fe/H] 



+0.05 (Hill et al. 2011). The resulting RC parameters are {{V — I)RCfl, ^/,_Rc)Buige = 
(1.09, —0.12). This derivation implicitly assumes that other parameters that effect HB 
morphology, such as age, a-enhancement, helium, and binarity, also behave in a linear or 
nearly-linear manner in the metallicity interval bracketed by 47 Tuc and NGC 6791. The 
potential impact of these uncertainties is discussed in Section |4.4[ 



4.3. (y — I)rc,qi Mj rc from Stellar Models, a Hipparcos calibration, and the 
Bulge RC Spectroscopic Metallicity Distribution Function 



Girardi & Safaris (2001) and Salaris & Girardi (2002) constructed a grid of horizontal 



branch stellar models to predict the functional dependence of the properties of the RC on 
various stellar population parameters. For the Galactic bulge, they combined the mean a- 
enhancement and metaUicity of [McWilham k Rich] ( [1994] ) ([Fe/H]= -0.22, [a/Fe]= +0.35) 
and assumed an age range of 8 to 12 Gyr to predict that the bulge RC should be 0.06 mag 
redder in {V — I) and 0.01 mag fainter in / than the stars of the solar neighborhood, which 



could be calibrated by using distances from the Hipparcos catalog (Ferryman et al. 1997). 



We update this estimate by using the same grid of stellar model outputs, and by as- 



suming the same age range. We use the metallicity distribution function of Hill et al. 



(2011), who measured 219 spectroscopic abundances in a sample of RC+RG stars with 



mean (V — /, /) corresponding to the measured position of the RC. They obtained a mean 
metallicity of [Fe/H]= +0.05, and mean [Mg/Fe]= +0.15, which we use as a proxy for the 
mean a-abundance. The 0.27 dex higher [Fe/H] and approximately half-as-high [a/Fe] shifts 
the bulge corrections, to 0.20 mag redder in {V — I) and 0.10 mag fainter in / relative to 
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that of the solar neighborhood. For the Hipparcos stars, the value of {V — I)rc,o ~ 1-01 
dPaczynski fc Stanek| [l998| ), and Mi^rc = -0.22 ± 0.03 dGroenewegeT^ [2058| ), yielding 
{{V — I)rc,o, Mj^Rc) = (1-21, —0.12). It is not clear why this result from theory is sub- 
stantially redder than our empirical estimates. It could be due to a difference in the mean 
metallicity of sightlines studied by Bensby et al. ( |2011[ ) and that studied by Hill et al. (2011 ). 
Another possibility is a different mapping between [Fe/H], the different a elements, age, and 
helium for bulge stars than that assumed by Girardi & Salaris (2001). 



4.4. Theoretical Estimates of the Effect of Population Uncertainties 

on Red Clump Parameters 

It is worthwhile to ask how the derived luminosity parameters of the bulge RC would 
vary if our population assumptions were to change. For this task we use the BaSTI stellar 



database ( Pietrinferni et al. 2004, 2006), whose breadth allows us to inspect the effect of 
several stellar parameters. The predicted effects of changing metallicity, age, and variable 
mass loss are computed using the synthetic HB generator, which interpolates between canon- 
ical tracks. We note that helium-enrichment is coupled to [Fe/H] as their models assume 
AY/AZ = 1.4. For helium-enrichment beyond (or below) this quantity, we must use the 
much coarser grid of He-enhanced models. We read off the difference in mean position of the 
HB tracks for the [Fe/H]= -0.70, t= 10 Gyr, Y= 0.30 and the [Fe/H]= -0.70, t=10 Gyr, 
Y=0.35 models. 

We thus obtain the following predicted differential relations for the brightness and color 
of the RC: 

Mi^RC oc +0.20[Fe/H] + 0.015(t/Gyr) - 0.029(100Y - lOOYss) + O.O6((5AM/O.1M0), (9) 

{V - I)RCfi oc 0.27[Fe/H] - 0.013(100Y - lOOYs^), (10) 

where Kjs is the scaled-solar helium abundance at metallicity [Fe/H], SAM parameterizes 
any potential difference in the mass loss relative to that assumed by the BaSTI database. 



Some readers may wonder why Equation 10 has fewer terms than Equation [9] for sufficiently 
high-metallicity, the mean color of the RC star becomes nearly independent of small changes 
in the total mass of the star, and thus either the mass-loss or the age of the progenitor. 



4.5. cri^RCfi from Stellar Models and the Spectroscopic Metallicity Distribution 

Function 



The intrinsic magnitude dispersion of the RC, cr j^rco^ is the magnitude dispersion that 
the bulge RC would have if the Galactic bulge were geometrically thin and if there was no 
differential extinction, i.e.: 



'^LRCfi — ^I,RC 



— (J, 



a 



where is the dispersion in distance modulus, and is the dispersion in extinction. 
Understanding of ai^RCfl is essential if one is to constrain cr^, a fundamental probe of Galactic 
structure ( Stanek et al.||1997 Cabrera-Lavers etaL][2008| . |Nataf fc Udaiskl] ( |20ll| ) estimated 
ci,Rc,o ~ 0.17, but that was done before realizing one of the findings of this work: that 
differential extinction can exceed values of 0.10 mag on scales as small as ~ 6' x 6' (see 



Section 6.1). We list three estimates: 



Girardi & Salaris (2001) used a large, detailed grid of stellar tracks, the metallic- 



ity distribution function of Mc William & Rich (1994), and a constant, metallicity 



independent star-formation rate over the range 8 < (t/Gyr) < 12 to estimate aj^RCfl 
0.107 mag. 



We use the synthetic HB calculator of the BaSTI stellar database ( Pietrinferni et al. 



2004, 2006) to estimate aj^RCfi = 0.031 for a simple, old, metal-rich stellar population. 



We quadratically add the effect of the dispersion in metallicity of the bulge (0.40 dex. 



Hill et al. 201 1|) and the predicted evolutionary effect from the BaSTl database of 

0.20 mag dex" 



variable metallicity on the RC luminosity, dMj/d[Fe/B. 
(^i,RC,o = 0.086 mag. 



to predict 



From Section 4.2, we have that (Jj^RCfi = 0.055 mag for a simple, old, metal-rich 
stellar population, and dMj/d[Fe/B.] = 0.11 mag dex~^. This yields an estimate of 
(^i,RC,o = 0.071 mag. 



The mean of the three estimates, dj^RCfl = 0.09 mag, is assumed in this work. 



4.6. (J — Ks)rc,o = 0.66 from Red Giant Color-Color Relations and Relative 

Reddening Expectations 



An estimate of {J — Ks)RCfi can be derived by combining the estimate (^ — /)ijc,o = 1-06 
with the numerous, precision measurements of RG color of Bessell & Brett (1988). We first 
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use the approximation that: 

K, = K + 0.1x{H ~ K). (12) 
we then interpolate rows 5 and 6 of their Table III to derive that: 

(J - K,)o = 0.620 + 0.625X ((F - /)o - 1.00), (13) 

yielding (J - Ks)RCfl = 0.66. 

A second estimate can be derived by studying the distribution of inferred E{J — 
Ks)/E{y — I) values (derived in Section |8| as the assumed values for (J — Ks)RCfl and 
{y — I)rc,q are varied systematically. For small errors in the intrinsic colors, denoted 

A(J — i^s)o = {J — Ks)o^lnierreA — {J — Ks)Q^True and A(\/ — /)o = ( V — /)o,Inferred — (^ — -?^)o,True: 

E{J-K,) EiJ-Ks) ^ r ^ AiV-I)o A(J-ir,)o \ ^^^^ 



E{V - I) ,nfe.re, E{V ' I) True V E{V ' I) E{J - K 

from which it follows that incorrect assumptions as to the values of (J — Ks)RCfi and {V — 
I)rc,o will cause the inferred values of E{J — Ks)/E(y — I) to be correlated with the 
reddening. We do a grid search and find that values of the intrinsic colors that remove the 
correlation between E(V^- J) and E{J - K,,)/ E{V - 1), i.e., cor{E{V -I),E{J -Ks)/E{V - 
I)) =0, satisfy the following relation: 

(J - Ks)o = 0.641 + OMOx{{V - /)o - 1.00). (15) 



For {V — I)RCfl = 1-06, Equation (15) predicts {J — Ks)rc,o = 0.66. Alternatively, requiring 
cor{E{V-I), E{J-Ks)/E{V-I)) = ±0.1 would shift iJ-Ks)Rc,o by 0.01, or iV-I)Rc,o by 
0.03, where a lower value of (J — Ks)rc,q or a higher value of {V — I)rc,q would necessitate a 
lower value of E{J — K^) / E{V — I) at higher E{V — I). We show later in the text that larger 
absolute values of the correlation are unphysical, and thus this method should determine 
(J — Ks)RCfl = 0.66 to an accuracy of 0.01, for a given value of {V — I)rc,o- 

A similar calculation can be performed to ascertain which values of (J — Ks)rc,o and 
{V — I)RCfl minimize the dispersion in E{J — Kg)/ E{y — I). The dispersion is minimized 
along the relation: 

(J - Ks)(i = 0.638 + 0.345X ((V - /)o - 1.00), (16) 

again yielding (J — -ft's)_RC,o = 0.66 for {V — I)RCfl = 1-06. However, the underlying physics 
of this estimate are very similar to those of the above estimate, so it is not an independent 
estimate, merely an enlightening reformulation of the second estimate. 



-17- 



5. The Reddening Law in V, I, J and Kg-. Theoretical Expectations 



Cardelli et al. (1989) combined data from several sources and found that the extinction 



law over the wavelength range 3.5/im > A > 0.125 /xm could be parameterized by a single 
variable, Ry = Ay/E{B — V). Thus, if one knows Ry and a reddening value toward a 
particular sightline, one can derive the extinction for each wavelength in that calibrated 
range. The equations of Cardelli et al. (1989) were updated by O'Donnell (1994), who 
obtained additional data. 



We convolve V and / filters of the Landolt photometric system (Landolt 1992) that 



OGLE-llI photometry is calibrated on (Szymahski et al. 2011) with a 4700 K blackbody 



curve, typical of bulge RC stars (Hill et al. 2011 ). We obtain effective wavelengths of 0.546 /xm 



and 0.804 fj,m for the V and / filters. For the J and Kg filters, we do as Gonzalez et al. (2012 ) 



and Indebetouw et al. (2005) , and respectively adopt 1.240 fim and 2.164 /xm. These effective 



wavelengths are a little different than those adopted by Schlegel et al. ( 1998 ), as those authors 



assumed the spectral energy distribution of an average elliptical galaxy for the source. 



The standard value of Ry for the interstellar medium is Ry = 3.1 (Schultz & Wiemer 



1975 Sneden et al. 1978). With the parameterization of Cardelli et al. ( |1989 ), the predicted 
reddening terms are Rj = Ai/E{V-I) = 1.424 and Rjkvi = E{J - K^) / E{V - 1) = 0.407, 
whereas they are Ri = 1.481 and Rjkvi = 0.416 with the parameterization of O'Donnell 
(1994). Given these numbers, we adopt Rj = 1.45 and Rjkvi = 0.41 as the standard values 



to which we compare our results. 

We comment on two additional uncertainties. The first is that due to instrumental 
uncertainty in the effective central wavelengths of the filters: a change of 1 nm in the effective 
wavelengths, approximately equivalent to a temperature change of 750 K, would modify the 
value of Rj by ~0.01 and Rjkvi by ~0.004. The second is that the effective wavelengths of 
the filters are modified when convolved with a foreground extinction. To gauge this effect, we 
assume Ry = 3.1, the equations of O'Donnell (1994), and apply 5 magnitudes of extinction 
in V. We find that Rj drops by ~0.02 and Rjkvi by ~0.007. These effects are real, but not 
large. 



6. Reddening Derived from the Mean Color of the Red Clump 

The reddening can be derived by taking the difference between the observed and intrinsic 
value of {V — I)rc- 



E{V -I) = {V- I)nc -{V- I)rc,o = {V - I) 



RC 



1.06. 



(17) 
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The distribution of reddening for 9,014 sighthnes is shown as a color-coded map in the top 
panel of Figure [5] 



6.1. The Scale of Differential Reddening 



Differential reddening, though less investigated than mean reddening, is of great impor- 
tance to Galactic bulge studies. Spectroscopic investigations of red giants assume photomet- 



ric temperatures and gravities (Fulbright et al. 2006; Zoccali et al. 2008 De Propris et al. 



2011 Hill et al. 2011 Johnson et al. 2011), and thus knowledge of the scale of differential 



reddening is a prerequisite to making more detailed analyses of the bulge giant metallic- 
ity distribution function. Microlensing investigations also rely on colors to constrain the 



properties of planet-hosting stars (Albrow et al. 1999 Yoo et al. 2004). 



For each of our fields, we have estimated the color dispersion of the RG-t-RC at the 
luminosity of the RC, (T[v-i),rc- This quantity will be equal to a quadratic sum of an 
intrinsic color dispersion and a reddening dispersion: 



{V-I),RC 



I z 

{V~I),RCfl "I" '-'^E(V-I) 



;i8) 



The distribution of color-dispersions looks remarkably similar to the distribution of 
E{y — I) shown in the top panel of Figure Isj We thus assume a simple parametric form for 



We plot the distribution of in the bottom panel of Figure [5 



a, 



{V-I),RC- 



(J, 



{V-I),RC 



^{V-I),RCfl + 



^2 
^DRl 



+ c 



DR2 



0.01 deg^ 



xE{V -If 



(19) 



where Afl is the solid angle subtended by the sightline, Cdri is the proportionality constant 
between total and differential reddening that is independent of the angular size of the sight- 
line, and Cr)ji2 parameterizes the component of differential reddening that depends on the 
solid angle of the field. We do a three-parameter optimization and find <J(v-i),rc,o = 0.121, 
Cdri = 0.070, and Cdr2 = 0.017, where the fitting was restricted to the 94% of sightlines 
with cr(^Y_j) jic < 0.25. As the mean solid angle for our sightlines is 0.011 deg^, our results 
demonstrate that significant differential reddening occurs on scales much smaller than a few 
arcminutes - it is essentially granular when reddening is measured with the RC method. 
The sum of Cdri and Cdr2 shows that to first order differential reddening averaging ~9% 
of total reddening, even over the small sightlines used by our investigation. 
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Fig. 5. — TOP: Color-coded reddening map of the Galactic bulge as observed by the OGLE- 
III VI photometric survey. BOTTOM: Distribution of o'(v-i),rc as a function of direction. 
Color-coded map of the color dispersion, cr{v-i),RC, which is the quadratic sum of the in- 
trinsic color dispersion and the differential reddening. The distribution of <7(v-i),rc looks 
remarkably similar to the distribution of E{V — I) shown in the top panel, suggesting a 
functional dependence of differential reddening on total reddening. For both panels, each 
color codes approximately equal areas. 
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6.2. Distribution of Dust Toward the Inner Galaxy 



We estimate the distribution of dust toward the Galactic bulge assuming a simple 2- 
parameter model: the mean density of dust along the plane pr, and the scale height Ho. 
Thus, the prediction for the reddening toward a given direction is: 



E{V-I){l,h) 



p^exp 



-rsm{\b\)/HD 



dr, 



where: 



Ro sin(a;) 



(20) 



(21) 



cos(6) sin(/ + a) ' 

We assume Ro = 8.4 kpc dGhez et al.|2008D and a = 25° ( |Rattenbury et al.|2007a| ), consistent 



with values estimated later in this work. 

We minimize the sum of the squares of the differences between the predicted and mea- 
sured reddening. This approximation yields pd = 0.427 magkpc"^ and Hn = 164 pc. The 
resulting scatter is 0.25 mag. This value of the scale height is larger than the 125 pc reported 



by Marshall et al. (2006). This is likely due to the fact a two-parameter model can only go 



so far. For example, pd may be a function of Galactocentric radius. 



6.3. Comparison to the All-Sky Reddening Map of Schlegel et al. (1998) 



Schlegel et al. (1998) produced a full-sky reddening map that is now one of the most 



widely used tools in astronomy. It is thus an important benchmark for comparisons. 

For each of our E{V — /)measurements we took the ratio to the nearest 4 pixel interpo- 



lation of E{B — V) from Schlegel et al. (1998), which we denote Rvi,BVsfd- We find a mean 
value of Rvi BVsid = 1-179 with a standard deviation of 0.173, the distribution is shown as 



a color-coded map in the top panel of Figure [6j The equations of O'Donnell (1994) predict 
that Rvi,BVsid = 0.24:9Rv + 0.562 , where we have made sure to use the filter definitions of 



Schlegel et al. ( 1998 ) for calculating E{B — V), while maintaining our definitions for E{V—I). 
The mean value of Rvi,BVsid = 1-179 implies Ry = 2.47, which is consistent with the value 
of Rv ~ 2.5 inferred later in this work. In the mean, the E{B — V) values of Schlegel et al. 
(1998) work well toward the bulge as long as one does not assume an Ry = 3.1 extinction 
curve. In the latter case, the predicted value of Rvi,BVsid = 1-33 would yield a systematic, 
~13% overestimate of the reddening. 

The 15% scatter in Rvi,BVsid does not seem to correspond to any known reddening index. 
Its correlation with E{J — Ks)/E{V — J) is p = —0.04, and its correlation with E{V — I) 
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Fig. 6. — TOP: Distribution of Rvi,bvsIA = EiV — I)/E{B — V)sid as a function of direction. 
The mean value of Rvi,BVsid = 1-179 is consistent with a mean extinction curve of Ry = 2.47, 
but the scatter of 0.17 does not correlate with any known reddening index. BOTTOM: Color- 
plot of distribution of measured dAj /dE{V — I) values as a function of position. Each color 
represents an equal surface area on the sky. More standard values of the reddening law, 
dAj/dE{V — /) > 1.35, are observed further from the minor axis. 
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is p = +0.02, both negligible values which undermine the prospect for an extinction-based 
explanation to the dispersion in -Ry/.BVsfd- 



7. The Reddening Law Toward the Bulge is Non-Standard and Non-Uniform 



There are two broadly used methods to convert a measurement of reddening into a 
derived extinction. The first is to assume a universal total-to-selective extinction ratio, i.e.: 



Ai 



E{V - I) 



xE{V-I) = RixE{V - I), 



(22) 



where we showed in Section |5] that Rj ~ 1.45 if one assumes standard extinction. 



The second method is to compute the linear regression of magnitude as a function of 
color, and to infer the total-to-selective extinction ratio in that manner: 



Rt 



A, 



dAj 



(23) 



E{v-i) dE{y-iy 

This method has the desirable characteristic that it is independent of any assumption of 
the intrinsic luminosity of whichever standard candle is being used: it depends only on 
the differential colors and magnitudes between sightlines, a robust quantity. This method 



(Udalski 


2003b 


Sumi 


2004 



Revnivtsev et al. 2010 Pietrukowicz et al. 2012), the near-IR reddening toward the bulge 



(Nishiyama et al. 2006, 2008, 2009), and to calibrate the (near-l-mid)-IR reddening in the disk 
dlndebetouw et al.||2005[ [Zasowski et"aL]|2009| |Gao et al.|[2009j ), among other studies. The 



assumption of Equation (23) is that the functional dependence of the variation in extinction 
on the variation in reddening is equal to the total-to-selective extinction ratio. The ~100 
deg^ Galactic bulge component of the OGLE-III survey allows us to ascertain whether or 
not this assumption is valid. 

We measure dAj/dEiy — /) as a function of position as follows. We keep the 9,014 
RC measurements qualified as reliable in Section [3| but remove those toward double-RC 
sightlines: (|/ + 0.5| < 4.5°) and (|6| > 4.75°). We then measure, for the 1,690 OGLE-III 
subfields used in this work that are not toward a double-RC sightline, the linear regression 
of {Irc + 0.03(/ — /central)) vs {V — I) RC of cvcry RC centroid measurement within 30' of 
the subfield central coordinate. The factor of 0.03(/ — /central) is a first-order correction for 
the Galactic bar's orientation which negligibly impacts the final values. The distribution 
of dAi/dE{y — /) as a function of direction is plotted in the bottom panel of Figure [6] 
Several results are immediately evident. We confirm the finding that dAi/dE(y — /) is 
smaller toward the inner Galaxy, as we obtain a mean value of dAj/dEiV — /) = 1.215. 
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Also evident, however, is that for sighthnes further separated from the minor axis we see 
higher, more standard values of dAj/dE{V — /), a significant number of sightlines with 
/ > 5° are color-coded black, meaning dAi/dE{V — /) > 1.35. No Galactic bulge fields with 
values of dAi/dE{V — I) so close to standard were found in the OGLE-II investigations of 
Udalski] ( |2003b| ) and |Simi| ( |2004| , as these were all closer to the Galactic plane and/or to 
the Galactic minor axis, and thus more heavily probed the extinction properties of the inner 



Galaxy. These values were also not found by the investigation of Pietrukowicz et al. (2012) 



as OGLE-III had lower cadence toward sightlines further from the Galactic minor axis, 
preventing an analysis of RR Lyrae stars. The trend toward more "standard" reddening 
properties for sightlines further from the Galactic center may be linked to the finding of 



Zasowski et al. (2009), who found that the functional dependence of Galactic extinction 



in the mid-IR depends on Galactocentric radius, with shallower (greyer) values measured 
toward the inner Galaxy. 

The histogram of dAi/dE(V — /) is shown in Figure [7| We plot the distribution of 
dAj/dE{V — I) as a function of mean {V — I)jic in Figure [s] The fact that the reddening law 
is smaller than the local value of Rj ~ 1.45 is itself independent of the reddening. However, 
there is a small trend between the two variables: we also measure a small correlation between 
dAi/dE{V - I) and {V - I)rc oi p = -0.11. Four of the fits to dAi/dE{V - I) are shown 
in Figure |9| The variability of the reddening law is a robust result. 
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Fig. 7. — Distribution of measured dAj/dE{V — I) values. The mean is 1.215 and the 
standard deviation is 0.18, with a mean error of 0.09. The observed range of optical reddening 
laws toward the bulge, 1.0 < dAj/dE(y — /) < 1.5, is smaller than the standard reddening 
law of dAj/dE{V — /) = 1.45, and variable with spatial direction. 
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Fig. 8. — Scatter plot of dAi/dE{V — I) versus the mean {V — I)rc of the RC centroids 
from which dAi/dE{V — I) is measured. Thick black X's denote the medians to 50 bins. 





Fig. 9. — We show the scatter plots of Irc vs {V — I)rc for four of the directions measured 
in this section, as well as the derived value of dAi/dEiV — I) with error. Plots are organized 
such that dAi/dE{V — I) appears in ascending order, clockwise from TOP-LEFT. 
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7.1. Attempt at Deriving the Extinction From a Single Reddening Value 



We attempt to derive the extinction to the Galactic bulge using three different methods 
that might appear valid. The three methods are to assume that: 

1. i?/ = 1.45 everywhere, the "standard" value of Ri when evaluating the extinction fits 
of ICardeUi et all (|l989| and lO'DonneUI (|T994|) at the value Ry 



3.1; 



2. Ri = 1.215 everywhere (the mean value found in the previous section); 



3. Ai/E{V - I) = dAi/dEiy - J), as per Equation (23) 



We plot the resultant values of Irca = Mj^rc + as a function of longitude on the left 



panels of Figure 10 On the right panels we plot the residuals of Irco relative to a moving 



fit as function of dAj/dE{V — I)- It is clear from Figure 10 that each of the three methods 
predicts large variations in the dereddened magnitude of the RC toward the same longitude 
in the Galaxy, which we will subsequently demonstrate to be a sign of failure. 



The failure of the first method, evident in the top panels of Figure [10} is the least 
surprising. It not only fails, it fails spectacularly. Huge structures are present, stretching 
downwards from the bar at every longitude, contributing to a scatter of 0.110 mag. At 
/ = 0°, the best-fit value of Irca is 14.131 mag. If one assumes Mj^rc = —0.12 (Section 
4.2), the implied distance to the Galactic center is ~7.0 kpc, very much on the low side. 
The correlation coefficient between dAi/dE(y — I) and Irco is p = —0.32, demonstrating 
that there is information in dAi/dE{V — I) not used by this method that could improve the 
fit. Though the non-standard VI extinction toward the bulge is already firmly established 
in the literature ( |Gould et all2001t |Udalski|[2003bl |Sumi||2004[ [Pietrukowicz et al.||2012| , we 
further demonstrate it here as it is a critical point worthy of reiteration. 

The second method does much better. The scatter is the lowest of the three methods, 
at 0.070 mag per point. The implied distance to the Galactic center is a reasonable 8.19 
kpc. However, the method still has some problems. The correlation coefficient between 
dAi/dE{V — I) and Irco is p = —0.23, demonstrating that there is pertinent information in 
dAi/dE(y — I) that could be used to improve the fit. There are various unphysical features 
in the plot, including a streak of points extending downwards near / = 0°, and two streaks 
extending upwards at / = +3°, +6°. 

Perhaps the most surprising result of this section is that the assumption that Aj/E{V — 
I) = dAj/dE{V — I), shown in the bottom two panels, does in fact not work. There are huge 
features that appear as an unphysical zigzag pattern superimposed on the Galactic bar. The 
correlation coefficient between dAj/dE{V — /) and Irco is p = +0.78. Its absolute value 
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is more than twice as large as that from the other two methods. However, the correlation 
coefficient has the opposite sign, further conffiming that there is information in dAi/dE{V — 
I) even if it is not yet clear what that information is. 

That Ai/E{V — I) ^ dAi/dE(y — /) is a shocking result. These systematics prove 
that the functional dependence of the variation in extinction on the variation in reddening 
is not equal to the total-to-selective extinction ratio. This assumption is often used in 
the literature without even stating that it is an assumption, demonstrating the extent to 
which this assumption appears natural. We understand that many readers will be skeptical 
about this point - it is not unreasonable to look at the line fits in Figure [9} to see how 
well they go through the points, and to conclude "that has to be extinction law". We 
dedicate the following subsection to the requirement of providing a heuristic explanation to 
this phenomenon. 



7.2. Variations in the Reddening Law Along the Line of Sight and Composite 

Extinction Bias 

The extinction law is a variable function of direction. Given that, we expect that it 
should vary within sightlines, and that the total extinction along a sightline to a distance R 
will therefore be an integral: 

where (dAj j dE(y — iy){r) and [dEiy — j dr){r) are the reddening law and reddening 



density as a function of radius. It logically follows from Equation [24] that the true extinction 
law for a given sightline will be an average weighted by the total reddening contributed from 
each type of interstellar medium and corresponding dust properties intersecting the line of 
sight. 

dAijdEiy — T) has historically been estimated by measuring the functional dependence 
of the variation in extinction on the variation in reddening. However, given that there will 
be many different extinction laws along a line of sight, that need not be the case. Consider 
a two-point regression, measured from RC centroids that are 30' apart on the sky. At a 
distance of 1 kpc, the two sightlines will be ~8 pc apart. At a distance of 4 kpc, the two 
sightlines will be ~30 pc apart. The difference in extinction will therefore be statistically 
biased toward whichever kind of dust is to be found deeper into the inner Galaxy, as distinct 
sightlines diverge linearly with distance. dAijdEiy — /) is a weighted average of different 
extinction laws, but the weights are not the same weights as those which go into Equation 
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Fig. 10. — Left panels show the scatter of dereddened RC magnitude, Irc,o = Mjjic + 
as a function of longitude for three different assumptions of the reddening law. Green 
line shows a moving linear fit to the points. Right panels show the residuals of Irc,Oi 
^iRCfi = -^i?c,o,Fit — Irc,o as a function of dAj/dE{V — I). Error bars denote the dispersion 
of AIrco ill each bin. 
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24 The name we assign to this phenomenon is "composite extinction bias" : The functional 
dependence of the variation in extinction on the variation in reddening is biased toward 
the kind of extinction which contributes differentially to different sightlines. A class of 
extinction that contributes equally to both sightlines will obviously contribute to the true 
value of Ai/E{V — I) of both sightlines, but it will not contribute to the scatter in Aj and 
E{V - /), and thus not to dAj/dE{V - I). 
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(l,b) = (0.5,-1.8) 
(14.67±0.02) + (1.00±0.01)xE(V-I) 

(l,b) = (0.5,-3.3) 
(14.43±0.02) + (1.24±0.01)xE(V-I) 




1 1.5 
E(V-I) 



2.5 



Fig. 11. — Scatter of Irc vs (V — I)rc near two sightlines. The measured reddening laws 
dAj/dE(y — I) are clearly different, but extrapolating the reddening laws back to E{V — I) = 
leads to an unphysical difference in the value of Irca + /i = 0.24 mag. We thus conclude 
that the slope dAi/dE{V — /) will be different for lower values of reddening. 



We show the derivations of dAi/dE{V — /) for two different directions in Figure 11 



These two fits yield different values of dAj/dE{V — /), and values of Irc,o that differ by 
0.24±0.03 mag. These two directions are at the same longitude, so the Galactic bar's ori- 
entation should affect both directions equally, and thus the difference in Irc,o = + Mj^rc 
cannot be due to a difference in distance modulus. It also cannot plausibly be due to a differ- 
ence in absolute magnitude, as RC stars have a metaUicity-dependence of ~0.2 mag dex~^ in 
/ (Girardi & Safaris 2001 Pietrzyhski et al. 2010). The difference of 0.24±0.03 mag is thus 
unphysical, but it directly follows from extrapolating the fits in Figure [TT] to E{V — /) = 0. 
We thus conclude that the fits should not be extrapolated back to E{V — I) = 0. There is a 
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component of the extinction that contributes to the extinction of these sighthnes, but much 
less strongly to the scatter of their extinctions. 

We are thus left with a dilemma. We have demonstrated that dAj/dE{V — /) 7^ 



Aj/E{V — I), yet the correlations found in Section 7.1 demonstrate that it has some physical 
meaning. There is no obvious way to probe any screen of extinction that would contribute 
equally to neighboring sightlines without making use of problematic assumptions as to the 
intrinsic values of distance modulus. We conclude that the conversion from reddening to ex- 
tinction cannot be done accurately with a single reddening value. The solution we propose 
is to utilize multiple reddening values, in different bandpasses. 



8. Combining E{V — I) and E{J — Kg) to Yield a Robust Extinction Estimate 
8.1. The Distribution of E{J - K,) /E{V - I) 



E{J 



For each of our sightlines with a measured RC centroid we add a measurement of 
Kg), in the manner described in Section 3.3 We thus have two completely inde- 
pendent reddening measurements for every sightline: E{V — /) and E{J — Kg). We plot 
the distribution of Rjkvi = E{J — Ks)/E{V — /) as a function of direction in Figure 12 



The mean value of Rjkvi is 0.3433, compared to Rjkvi ~ 0.41 if one assumes standard 
dust properties, as shown in Section |5] Once again, the dust toward the inner Galaxy is 
characterized by a steeper extinction curve than standard. Further from the minor axis, a 



reassuring similarity is observed between Figure 12 and the bottom panel of Figure |6| Both 
have (relatively) more standard extinction further from the Galactic minor axis, toward 
approximately the same directions of / = —4°, +5°, 6 ~ —5°. 

We note two advantages of using Rjkvi rather than dAj/dE(y — /). The first is that 
the ratio of two reddenings can be taken toward both single-RC and double-RC sightlines. 



as double-RCs do not differ in color (Nataf et al. 2010 McWiUiam & Zoccah 2010). The 



second advantage is more fundamental: whereas the measurement of slopes dAj/dE{V — I) 
must be done over ~30' scales to be precise, independent values of Rjkvi can be measured 
for each sightline, giving us very sharp resolution on the extinction law. The fact each of 



9,014 measurements shown in Figure 12 is independent of the other measurements confirms 
that the variations in the extinction law are real: if the variations were due to noise, they 
would not cluster together. 
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Fig. 12. — Color map of the observed distribution of Rjkvi = E{J — Ks)/E(y — /). If 
the extinction curve toward the bulge were standard we would measure Rjkvi ~ 0.41, in 
contrast, even the demarcation between the 6th and 7th septiles, Rjkvi = 0.37, is of a 
steeper extinction curve. The mean value, Rjkvi = 0.3433, implies Ry ~ 2.5 (Cardelli et 



al. 


1989 


O'Donnell 


1994) 



8.2. Aj As a Function of E{V - I) and E{J - Ks 



The parameterizations of Cardelli et al. (1989) and O'Donnell (1994) assume that if 
one knows Ry, then one knows all the total-to-selective extinction ratios in the optical and 
near-IR. We do not know Ry as we do not have measurements in 5-band, however, knowing 
Rjkvi is equivalent to knowing Ry if one assumes the extinction law is a single-parameter 



function. Both Cardelli et al. (1989) and O'Donnell (1994) parameterize the extinction curve 
as Ax = Ay{a{X) + b{X)/Ry), allowing us to obtain linear relations linking Aj, E(y — I) 
and E{J - K, 



We invert the equations of Cardelli et al. (|1989|) and find: 
Ai 



E{V-I) 



0.1713 + 3.078X 



E{J - Ks) 
E{V -I) ' 



(25) 



which is equivalent to: 



Ai = 1.228xE{V - I) 



1 + 2.507x{RjKvi - 0M33) 



(26) 
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whereas if we use the values of O'Donnell (1994): 



Ai = 1.266xE{V - I) 



1 + 2.323x{RjKvi - 0.3433) 



(27) 



There is a simple interpretation to Equations (26) and (27). They are first-order expansions 



to the reddening law in the expansion variable Rjkvi, with smaller values of Rjkvi implying 
a steeper extinction curve. A steeper extinction curve is conventionally interpreted as being 



due to smaller dust grains (Draine 2003) 



Figure 13 shows the result of three different methods to fit for the extinction to the 
Galactic bulge. The top two panels show the result for Aj = 1.215xii^(V^ — /), with no 
information from E{J — Kg), as a comparison. The correlation between the residuals to 
the moving fit and the dereddened apparent magnitudes is p = —0.39, a large number 
demonstrating the failure of this method to use all the information available. The scatter is 
0.070 mag. 

In the middle panel, we plot the result given an extinction prescription that is the 
mean prediction of ICardelh eTaL] (p89| and [O'Donnelll ([l994l): Aj = 1.2A7xE{V - /)(1 + 



2A15x(Rjj^vi ~ 0.3433)). Surprisingly, this does less well than not using the information 
from E{J — Kg): the Icr scatter is increased to 0.078 mag. The correlation between the 
residuals to the moving fit and the dereddened apparent magnitudes is p = +0.59. That is 
nearly twice as large as that obtained when simply using Ai = 1.215, though of the opposite 
sign. 



In the bottom two panels, we plot the results for the extinction prescription: 

Aj = 1.217xE{V - I) l + 1.126x(i?j^yj- 0.3433) 
= 0.7465xE(y - /) + 1.3700xE(J - K^), 



{2t 



values optimized by minimizing x^- We show the equation in two algebraic formats for clarity. 
This method produces the smallest scatter, at only 0.060 mag. The bottom-left panel shows 
the fewest structures. Finally, the correlation between the residuals to the moving fit and 
the dereddened apparent magnitudes is p = +0.11. It is the closest to zero of the three 
methods. That it is not exactly equal to zero may be due to the fact that the error on Irc is 
a function of location on the sky, and is therefore evidence that there are third parameters 
to the extinction law, the first two parameters being E{V — I) and Rjkvi- These third 
parameters could be factors acting on only certain types of dust grains, or they could be 
phantom extinction parameters such as gradients in metallicity, age, or helium, which could 
easily masquerade as a variation in the reddening law. 
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A,= 1.215xE(V-I) 
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Fig. 13. — Left panels show the scatter of dereddened RC magnitude, Irc,o = Mj ^c + 1^ 
as a function of longitude for three different assumptions of the reddening law. Green 
line shows a moving linear fit to the points. Right panels show the residuals of Ircai 
^Irc,o = -^RCo.Fit — Irc,o, as a function of E{J — Ks)/E(y — I). Error bars denote the 
dispersion of AIrco in each bin. 



Nevertheless, whatever these third parameters are, they are small. There are almost 
no features in the bottom-left panel of Figure 13, unlike all other methods we previously 
have shown in this work. The scatter of 0.06 mag is very close to the noise floor: given our 
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average errors of aj^^c = 0.034 mag, cr{v-i),Rc = 0.01 mag, and (J(j-Ks),rc = 0.01 mag, we 
could do no better than a dispersion of 0.039 mag. This noise floor does not include the 
impact of metallicity gradients and latitude-dependent projection effects. The method of 
Equation (28) yields a scatter of 0.06 mag, compared to 0.10 mag when assuming a standard 
reddening law of Aj = lA5xE{V — /)). That is a factor 4 reduction in the variance once a 
conservative estimate of the noise floor is removed. Thus, to that high degree of accuracy, 
we have solved the longstanding observational challenge of the non-standard VI extinction 
toward the Galactic bulge. 



8.3. Caveats to the Reddening Law 

We discuss four caveats to the reddening law that can induce second-order effects on 
Galactic bulge studies at the ~0.05 mag level. 

The first is that there is evidence for a dependence to the reddening law on the spectral 
energy distribution of stars being investigated. The reddening law found in this investigation 



is a bit shallower than that reported by Pietrukowicz et al. (2012), who used RR Lyrae stars 



as standard candles. This is partly because OGLE-III had lower cadence toward sightlines 
further from the minor axis, where the reddening law is relatively more standard. However, 
even once that is accounted for there is still a small effect. We compare RRab stars to the 
closest of the clean RC centroids described in Section [7] that are also not toward a double-RC, 
and that are within 5' of an RRab star. We find: 

(V - I)rc -iy- 1) RRab = 0.552 - 0.084X {{V - I)rc - 2.22) - 0.83 x{\h\Rc- \h\RRa, - 0.021), 

(29) 

The latitude term captures the fact that an RRab star will be more reddened than the 
nearest RC centroid if it is closer to the plane. No significant trend is found for Irc — iRRab 
regardless of whether or not we fit for a longitudinal term, so the effect is purely in V- 
band. This temperature-dependence is several times larger than that predicted from theory 
in Section [5j The origin of the significance to the terms in Equation |29] is as such unresolved. 

Metallicity gradients could have a small impact. From Section|4| we expect that for a 0.2 
dex~^ increase in metallicity at fixed age, from [M/H]= to [M/H]= +0.20, Mj^rc should 
get fainter by 0.04 mag, and {V — I)rc redder by 0.06 mag, thus making (J — Ks)rc redder 
by 0.04 mag. A 0.2 dex increase in the metallicity would therefore cause "dereddened" 
RC stars to appear brighter by ~0.03 mag, as the effect of overestimating the extinction 
(due to redder intrinsic colors) would be larger than that of the dimming of the RC in 
/. The actual amplitude will be modified depending on how [Fe/H] correlates to [a/Fe], 
age, helium abundance, and binarity. The bulge vertical metallicity gradient measured at 
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large separations from the plane is ~0.6 dex kpc ^ (Zoccali et al. 2008 Johnson et al. 



2011), though evidence suggests this flattens out within 4° of the plane ( Ramirez et al .112000 



Gonzalez et al. 2011c; Rich et al. 2012). An additional source of uncertainty is that the 



mapping of other abundances onto [Fe/H] may depend on kinematics (Johnson et al. 2012), 
and the Galactic bulge age-helium-metallicity relationship is only loosely constrained at high 
metaUicities (Nataf & Gould 2012). Once a clear picture of the bulge metallicity gradient 
emerges, in both latitude and longitude, there may be a need for an additional iteration to 
reddening and extinction maps. 




Fig. 14. — The scatter in AIrco = -^i?,c,o,Fit — Irc,o, as a function of four parameters. There 
are no strong trends with E{V — /), E{J — Kg) or E{J — Ks)/E{V — I). However, there are 
clear, non-linear trends with separation from the plane. For each panel, error bars denote 
the dispersion of AIf>cfl in each bin. 



There appears to be a latitude term in the residuals to the extinction law when using 
Equation (28). We show the scatter in Irc,o as a function of four parameters in Figure 



14 Whatever trend there may be with E{V - I), E{J - Ks) or E{J - Ks)/E{V - I) is no 



greater than ~0.02 mag over the bulk of the data. However, there is a large, non-linear trend 
with separation from the plane, represented as the absolute value of the latitude. At large 
separations from the plane, /rc,o is ~0.05 mag brighter than the faintest Irca values seen, at 
around |6| = 3°. This can plausibly be due to projection effects; at large separations from the 
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plane, the distance to the maximum density of stars along a line of sight through a Galactic 



bar becomes systematically smaller than distance to the major-axis (Cabrera-Lavers et al. 



2007). If that is the explanation, the same offset will show up in studies of the RC in Ks and 



in mid-IR filters. However, for \h\ < 2°, we see another very large increase, where the RC 
becomes up to ~0.1 mag brighter than its mean longitudinal value. Further investigation of 
this feature is warranted. An additional parameter to the extinction law may seem the most 
likely culprit, but variations in the metallicity or projections effects could also contribute. 
An extension of our methodology to incorporate additional reddenings, such a,s E{B — V), 
E{I — J), or mid-IR reddenings could help disentangle whether the remaining scatter is due 
to additional parameters of the properties of the interstellar medium, or due to the intrinsic 
properties of the bulge stellar population. 
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Fig. 15. — The predicted relationship for Aj/E{V — /) vs E{J — Ks)/E{V — I) using the 
equations of O'DonneU] ( |1994 ) (blue) as written, using the equations of O'Donnell (1994) 
modified to assume A{X) oc A~^ °° to be consistent with the result of Nishiyama et al. 



2009 



(red), and the best-fit relationship found in this work (black), with binned medians shown as 
black circles. The red and black curves intersect at {E{J — Kg)/ EiV — /), Ai/ EiV — /)) = 
(0.337,1.208), very close to the respective mean values of 0.3433 and 1.217 found in this 
work. Histograms of E{J — Ks)/ EiV — I) and Ai/E{V — I) values overplotted. 
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Finally, we note what may be the most important caveat to the reddening law. Our 
mean reddening values, Ai/E{V — I) = 1.22 and Rjkvi = 0.3433, both suggest an Ry ~ 2.5 
extinction curve. However, this is done without any actual S-band measurements, and thus 



the inference could be affected by systematics. The formalism of Cardelli et al.| ([1989) and 
O'Donnell (1994) both assume a universal extinction law in the near-IR, A{X) oc A~^'^^, 
that is independent of Ry- However, it has been demonstrated that the near-IR reddening 



law toward the inner Galaxy is steeper. Nishiyama et al. (2009) used RC stars to measure 
A{X) oc A^^°°, a finding supported by Fritz et al. (2011), who used line-emission to study 
the extinction curve. 

It may be that the exponential slope of the near-IR extinction law rises to higher values 
for lower values of Ry- Cardelli et al. (1989) and O'Donnell (1994) would not have been 



able to measure this, as none of their 51 sightlines probed values of Ry as low as the mean 
value toward the bulge inferred in this work. This could explain why the coefficients to 



Equation (|28|) do not match the predicted values. In Figure 15 , we show that adjusting the 

-2.00 yigifig predictions more consistent 



equations of O'Donnell (1994) to assume A{X) oc A 
with our measurements for the typical values of E{J — Kg) / E{V — I) found in this work. This 
may be a coincidence, but it would be worth investigating, once data in more bandpasses 
become available, whether the exponential slope of the near-IR extinction curve rises as Ry 
goes down. We note that studies of the extinction in the hosts of type la SNe routinely 
find an Ry ^ 2.5 extinction curve (Guy et al. 2010 Chotard et al. 2011 Mandel et al. 



2011), and then subsequently assume the functional dependence on wavelength calibrated in 



investigations of Milky Way extinction. If that functional form were modified, it would have 
significant implications not just for Galactic studies, but also for cosmology, 



9. Why Studies of RR Lyrae in MACHO Photometry Have Not Constrained 

the Reddening Law Toward the Bulge 



Kunder et al. (2008) investigated the reddening law toward Galactic bulge RR Lyrae 
stars observed by the MACHO survey. They inferred a reddening law of RyyR = Ay/E{V — 
R) = 4.3 ± 0.2, consistent with the standard extinction curve resulting from Ry = 3.1. It 
is the only investigation in the literature arguing for a standard reddening law toward the 
Galactic bulge, and thus needs to be understood. We demonstrate here that this result is 



likely to be a consequence of an incorrect assumption: that the transformations of Alcock et 



al. (1999) correctly calibrate the MACHO V and R filters into Johnson V and Kron-Cousins 
R filters. 




Fig. 16. — Comparison of RR Lyrae stars from MACHO ( Kunder et al.|2008 ) with those from 
OGLE-III ( Soszyhski et ah [2011 ). TOP: The mean magnitude of RR Lyrae stars in Vogle-iii 
is only 0.16 mag brighter than the magnitude at minimum hght in ^acho? demonstrating 
that Vmacho is non-standard. BOTTOM: {V — /)ogle-iii — (V — -R)macho as a function of 
E(V-I). Extrapolating to E{V - I) = yields {V - /)ogle-iii - (V - R)macuo = 0.014, an 
impossible result if MACHO photometry were standardized. 



We test this by first matching the MACHO RR Lyrae catalog of Kunder et al. (2008) 
and OGLE-III RR Lyrae catalog of Soszyhski et al. (2011 ), keeping only the RRab stars, for a 



cross-match of 2000 sources. We plot the comparison in Figure [16} In the top-panel, we show 
that the mean value of Vogle-iii is only 0.16 mag brighter, on average, than the minimum 
light value of \4iacho- However, the mean \^-band amphtude of RR Lyrae stars reported 
by Kunder et al. (2008) is 0.88 mag. We have downloaded some RR Lyrae lightcurves 
from AS AS (Pojmanski 2002), and found that the magnitude difference between minimum 
light and (time-averaged) mean light is typically 35-50% of the amplitude of variation, for a 
difference of ~0.3 mag in V between minimum and mean brightness. The 0.16 mag difference 
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is therefore unphysical if both Vogle-iii and Vmacho are standard. We conclude that \4^acho 
is centered at a lower effective wavelength than standard l^-band, even after applying the 
corrections of Alcock et al. (1999). The offset is very large, equivalent to ~0.15 mag. 



In the bottom panel of Figure 16 we plot the behavior of (V^ — /)ogle-iii — (V^ — -R)macho 
as a function of E{V—I). If the filters were standard, this would effectively be a plot of (R—I) 
vs E(y — /), but the results of the fit demonstrate this is not the case. The slope obtained 
for (V^ — /)oGLE-iii — — -R)macho vs E{V — I) is 1.01, much too high to be due to the effect 
of the reddening law on standard filters, as Ar ~ 0.5x (AyH- A/). Moreover, the extrapolation 
of the line back to E{V — /) = yields {V — /)ogle-iii — {V — -R)macho = —0.02 mag. 
These dereddened colors are impossible for standard filters. They are inconsistent with the 
assumptions that {V — -R)macho,o ~ 0.28 (iKunder et al.||2008|) and {V — /)ogle-iii,o ~ 0.47 



(Pietrukowicz et al. 2012). 



The OGLE-III filters are rigorously calibrated on the Landolt photometric system to an 
accuracy of a few hundredths of a magnitude ( Szymahski et al.||2011 ). We therefore conclude 
that there is an error in the calibration of the MACHO filters. 



10. Constraining Fundamental Parameters of Galactic Structure I: The 
Galactocentric Distance and the Viewing Angle to the Galactic Bar 



The structure of the inner Galaxy is currently a matter of active research, as the spatial 
morphology of the bulge is not fully known. There is an X-shaped component at large 



separations from the plane (Nataf et al. 2010 Mc William & Zoccali 2010; Saito et al. 2011 



Li & Shen 2012 Ness et al. 2012), evidence for a long bar at large separations from the 



minor axis (Benjamin et al. 2005 Cabrera-Lavers et al. 2008) that may be due to leading 



ends ( Martinez- Valpuesta &: Gerhard|2011 ), and an inner structure that is either a secondary. 



nuclear bar or a viewing effect (Alard 2001 Nishiyama et al. 2005 Gonzalez et al. 2011b 



Gerhard & Martinez- Valpuesta 2012 Robin et al. 2012). There is also evidence that the 



metal-poor stars are distributed as a classical bulge ( Babusiaux et al.||2010 De Propris et al. 



2011 Uttenthaler et al. 2012), whereas the data on metal-rich M-giants are consistent with 



entirely cold kinematics (Kunder et al. 2012). The functional dependence of RC brightness 



on direction has historically been a powerful constraint on bulge properties, responsible for 
many of the insights mentioned here. It is thus interesting to see what this new dataset of 
dereddened RC centroids can tell us about the geometry of the bulge. 

We plot the dereddened apparent magnitudes to the RC centroids measured in this 



work in Figure 17, where we assume Mj^rc = —0.12, as per Section 4.2, We use the same 
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Galactic coordinate system as Dwek et al. (1995). A line fit in the X-Y plane to RC centroids 
satisfying |Z| < 3.0° yields a Galactocentric distance Rq = 8.20 kpc and an apparent viewing 
angle of a = 40°, which is a soft upper-bound and consistent with a true viewing angle of 
a = 25 — 27° (explained below). 




X (Kpc) 



Fig. 17. — Distances to dereddened RC centroids are projected onto a face-on view of the 
central region of the Milky Way. Black points show the projected location of each measured 
RC centroid, and green points show 50 binned values. 



Our estimate of Rq = 8.20 kpc resolves a dissonance in the literature. The extinction 
map of Sumi (2004), which assumed that Aj/EiV — I) = dAj/dEiV — I), yielded a distance 
to the Galactic center of ~9 kpc (Rattenbury et al. 2007a; VanhoUebeke et al. 2009) when 
applied to VI photometry of RC stars. That value is too large when contrasted to the 
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geometrically determined distances to the Galactic center of 7.94 ± 0.42 kpc (Eisenhauer et 



al. 2003), 8.27 ± 0.29 kpc (Schonrich 2012), and 8.4 ± 0.4 kpc (GhezetaL 2008). It is also 



larger than distances derived from studies of the RC in the near-IR (Babusiaux & Gilmore 



2005 Nishiyama et al. 2006). Rattenbury et al. (2007a) actually shifted their measured values 
of Irc by 0.3 mag to force the bulge to be centered at i?o = 8 kpc, whereas [Vanhollebeke et 



al. (2009) commented that "The origin of this discrepancy is not understood at the moment". 



Our distance of 8.2 rather than 8.0 kpc reducehs that offset by 0.05 mag, and our corrected 
calibration og Mj^nc further reduces it by 0.14 mag. The different extinction law accounts 
for the rest. The solution to the non-standard VI extinction toward the Galactic bulge thus 
relieves a major bottleneck in Galactic bulge studies: there is no peculiar population effect 
or catastrophic failure of stellar evolution models decalibrating our primary standard candle 
by a spectacular value of 0.3 mag. 

The margin for the population correction of the bulge RC, and thus the error on our 
estimate Rq = 8.20 kpc, is rapidly shrinking. The spectroscopic metallicity distribution 



function of the Galactic bulge is now known to an impressive degree of accuracy (Zoccali et 



aL] |2008| |Alves-Brito et al.|[20T0| [Hill et al.|[20TT| [Bensby et aT][20TT| [Johnson et aL][20TT 



Rich et al. 2012). The same is not true of the age and helium abundance of the Galactic 



bulge. As the uncertainty in these two parameters is only an issue for the more metal-rich 



stars (Clarkson et al. 2011 Bensby et al. 2011 Nataf & Gould 2012), its integrated effect 



will not be very large. A lower age or higher helium abundance for the bulge would both 
require a larger value of Rq. 

The value of a = 40° is a soft upper bound. This is because the distance along the 
plane to the maximum density along the line of sight to a triaxial structure is strictly less 
than the distance to that structure's major axis for that sightline on the far side, and strictly 
greater on the near side. The difference can lead to a bias of ~50% in the inferred viewing 



angle of the Galactic bar, and depends on the bar's axis ratios (Stanek et al. 1994 Cabrera- 



Lavers et al. 2007). This can also be discerned by investigating viewing effects by means 



of sophisticated dynamical models. We compare our results to predictions from the N-body 



model used by Martinez- Valpuesta & Gerhard (2011) for two latitudes in Figures 18 and 
19 As their analysis assumes R=8 kpc, we have rescaled their distances by 1.025 to be 
consistent with the distance inferred in this work. Though their assumed viewing angle is 
a = 25°, their model points line up nearly perfectly with our data points. The apparent 
angle from the model is ~ 1° greater for the data in Figure 18 and ~ 5.5° greater in Figure 
19 As the apparent angle is twice as large as the true angle, this suggests a superior match 
would be attained if the assumed angle of the model was set to a ~ 27°. 
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Fig. 18. — Data from this work (black), binned data from this work (green), and predictions 
(magenta) from the N-body model of Martinez- Valpuesta & Gerhard (2011), for sightlines 
satisfying (-10.5° < / < +5.5°, 2.5° < |6| < 3.5°). 




X (Kpc) 



Fig. 19. — Color scheme as in Figure 18 , for sightlines satisfying 
\b\ < 3.5°). 



-10.5° < / < +10.5°, 2.5° < 
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10.1. The Apparent Magnitude of the Red Clump as an Improved Constraint 

on Microlensing Events 

A fundamental parameter of microlensing lightcurves is the angular size of the Einstein 
ring, which is proportional to the square root of the mass of the lens. It is related to 
observables by: 

Be = -, (30) 
P 



where p is a term due to finite-source effects that can be directly measured from well-sampled 
lightcurves ( |Gould||1994t [Nemiroff fc Wickramasinghe||1994[ |Witt fc Mao||1994D . The angular 
size of the source, 9^, is given by: 

F = nSe,^, (31) 
where F is the dereddened flux of the source, and S is the dereddened surface brightness. 
The first step used in inferring the surface brightness is measuring the color difference 



between the source and the RC, where the intrinsic color of the RC is calibrated by Bensby 



et al. (2011). A combination of color-color relations (Bessell & Brett 1988) and empirically- 



calibrated relations between color and surface brightness (Kervella et al. 2004) can then be 



used to yield the surface brightness. It is not as straightforward to infer the flux, due to 
the lack of a calibration for the intrinsic magnitude of the RC and the degeneracy between 
extinction and distance. The flux has been inferred by assuming a value of Mjnc, a distance 
to the Galactic center, and a correction for the distance due to the Galactic bar's orientation 



angle (Yoo et al. 2004). 



The results of this work will greatly simplify this process by streamlining the step that 
may have been the largest source of error: the determination of the dereddened apparent 



magnitude of the source. The bottom-left panel of Figure 13 shows that the dereddened 
apparent magnitude of the RC, Irc,o, is a very well-behaved function of longitude, for which 
we estimated an intrinsic scatter of ~0.04 mag in Section 8.2 We summarize the results 



in Table [TJ The need to assume a distance to the Galactic center, a value of Mj^fic, and a 
correction due to the Galactic bar's orientation is no more: one can now simply measure the 
brightness difference between the RC and the source, and then read off the value Irca as a 
function of / from Table [H 
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Table 1: Dereddened magnitude of the RC, Irc,o, as a function of Galactic longitude. 
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11. Constraining Fundamental Parameters of Galactic Structure II: 
The Stellar Mass of the Galactic Bulge 

The stellar mass of the Galactic bulge is currently unknown. A widely quoted value is 
I.SxIO^^Mq from the work of Dwek et al. (1995), who converted the bolometric luminosity 



In contrast, 



of the bulge into a number of RG stars, and then obtained a total mass by integrating a 
Salpeter initial mass function (Salpeter 1955) down to M = O.IMq 
estimated 



1995 



1996 



'1.7 — 2.8)xlO^''M0 by using the tensor virial theorem, and Zhao & Mao 
estimated M > 2.OxlO^'^M0 based on the frequency of gravitational microlensing 



Blum 



events toward the bulge. Robin et al. (2012) estimated I.IxIO^^Mq by fitting 2MASS data 



to various parametric models. 

A baryonic mass for the Galactic bulge would constrain models of Milky Way assembly, 
both in its own right and by anchoring the mass scale of the Galaxy as a whole. It is thus 
interesting to ask how well this dataset could eventually contribute. Most RC studies have 
been focused on tracing distances, but the Paczynski-Stanek equation also allows one to 
trace star counts, as Nrc is a free-parameter. We thus provide a constraint on the stellar 
mass of the Galactic bulge by estimating the number of RC stars. 



We plot the surface density of RC stars, S/jC; in the top panel of Figure [20} where we 
normalize to the density of 45,100 RC stars deg~^ observed toward Baade's window. We also 
show the scatter of T,rc for latitudinal and longitudinal stripes on the top-left and top-right 



panels of Figure |2T| respectively. The top-left panel demonstrates that S^c" is maximized 
near / = 0° for all latitudinal stripes, however, the sharp peak turns to a plateau at large 
separations from the plane. The top-right panel shows the expected result that number 
counts increase monotonically with decreased separation from the plane at fixed longitude. 

We measure 2.94x10^ RC stars in OGLE-111 over a viewing area of 103.75 deg^ with 
an overestimate of ~10% in each due to double counting arising from overlapping fields. 



We extrapolate this across the bulge using the "Gl" model of Dwek et al. (1995), assuming 
a viewing angle a = 25°, a corotation radius of 4.0 kpc, and a galactocentric distance of 
Rq = 8.20 kpc. We find a best-fit predicted RC population for the bulge of 12.7x10^, with 
output axial ratios of Xq : Yq : Zq = 1 : 0.45 : 0.30. In contrast, fixing a = 40° yields 
12.3x10^ RC stars with axis ratios of Xq : Yq : Zq = 1 : 0.35 : 0.35. The latter fit is less 
optimal, by Ax^ = 211. 
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Fig. 20. — TOP: Surface density of RC stars toward the Galactic bulge Tijic, as a function of 
direction. Values are normalized to the surface density toward Baade's window (/ = l°,b = 
—3.9°) of 45,100 RC stars deg~^. BOTTOM: Distance modulus dispersion of bulge RC stars 
as a function of direction, after application of a 20' smoothing.. 
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Fig. 21. — The number density of RC stars (normalized to that of 45,100 deg~^ toward 
Baade's window) and dispersion in distance modulus (with 20' smoothing) as a function 
of direction. Panel (a) Tjrc as a function of longitude for \b\ ~ 2.5° (black), |6| ~ 4.0° 
(blue), and |6| ~ 5.5° (red). Panel (b) function of longitude for three latitudinal 

stripes, same color scheme as panel (a). Panel (c) function of latitude for / ~ 0.0° 

(magenta), / ~ 5.0° (cyan), and / ~ 10.0° (grey). Panel (d) T^jic as a function of latitude 
for three longitudinal stripes, same color scheme as panel (c). 



Converting this population number into a stellar mass requires assumptions as to the 
stellar lifetime of the RC and the IMF of the bulge stellar population. For an old, metal-rich 



population, Nataf et al. (2011c) used stellar models to compute that the initial mass of stars 



found on the RG branch (and thus at the helium flash) is predicted to be: 

log(M/M0) = 0.026 + 0.126X [M/H] - 0.276x log(t/(10 Gyr)) - 0.937x (Y - 0.27). 



(32) 



The lifetime of the helium-burning phase is approximately 100 million years (Salasnich et 



al. 2000). Assuming a mean metallicity of [M/H]= -1-0.16 (Hill et al. 2011), a mean age of 
t = 10 Gyr dClarkson et~aL][20lIj [Bensby et aL][20lI| ) and a mean initial helium abundance 



of Y=0.33 ( [Nataf et aL]|2011a| [Nataf fc Gould] |2012D leads to the estimate that RC stars 
originate from the initial mass range (O.9744M0 < M < 0.9771Mq). We assume a Salpeter 



IMF (Salpeter 1955) for the bulge over the range 0.1 Mq < Minitiai < WOMq, that 90% of 



bulge stars that complete their stellar evolution through the helium-burning phase (discussed 



below), and the same remnant mass function as Gould (2000), whose assumptions are based 



on the results of Bragaglia et al. (1995) and Thorsett & Chakrabarty (1999). We end up 



with an estimated stellar mass for the Galactic bulge of 2.OxlO^"M0. 

There are several uncertainties in this estimate, which are tractable in principle. The use 



of a triaxial model from Dwek et al. (1995) is problematic - these triaxial ellipsoid models are 



now ruled out by the data, for reasons listed at the top of Section 10 However, our purpose in 



this section is merely to provide an estimate. We publish all our measured RC parameters and 
it will be straightforward to reproduce our calculations once more sophisticated structural 
models are available. An additional underestimate arises from our assumption that 90% of 
stars that "should" have ended up on the RC did end up on the RC. In reality, many bulge 
stars either end up on the extreme blue horizontal branch or skip the helium-burning phase 
altogether. The HST photometry used by Clarkson et al. (2011) shows that at least ~3% of 



bulge horizontal branch stars are not RC stars, and the fraction could be found to be higher 
if deeper photometry is obtained. Moreover, 10% of field white dwarfs have masses lower 
than the helium-ignition limit ( Liebert et al.||2005 ), some of which are not in binaries (Brown 



et al. 2011). These would not even show up on the blue horizontal branch, but should be 



incorporated into a complete model. Our estimate assumes that the fraction of bulge stars 
that skip the helium-burning phase can be inferred from the fraction of local white dwarfs 
with masses smaller than that of the helium-burning limit. The third source of error is in our 



evolutionary assumptions for Equation 32 a higher stellar mass for the bulge would result 
from assuming a lower mean age, a lower helium abundance, or a higher mean metallicity. 



Finally, the assumption of a Salpeter IMF is likely the most significant source of error: 
our estimate of the total stellar mass would drop by one half if we assumed the same IMF 
as Robin et al. (2012). Zoccali et al. (2000) used HST observations of the bulge luminosity 
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function to estimate a = —1.33 ± 0.07, which is very bottom-hght relative to the value of 
a = —2.35 for a Salpeter IMF. Calchi Novati et al. (2008) used the duration distribution of 
Galactic gravitational microlensing events toward the bulge to estimate a = —1.70 ±0.5. On 
the other hand, Conroy & van Dokkum (2012) have argued that the properties of M-dwarf 
spectral features in the integrated light of field galaxies implies that stellar systems which are 
metal-rich and a-enhanced have bottom-heavy IMFs. Their analysis is not consistent with 
that of Zoccali et al. (2000) and Calchi Novati et al. (2008), and thus further investigation 
is warranted. The spectacular Galactic bulge luminosity functions of Brown et al. (2010), 
measured with HST, may provide a path to determining the Galactic bulge IMF, and thus 
the total bulge stellar mass. 



12. Constraining Fundamental Parameters of Galactic Structure III: 
The Geometrical Thickness of the Galactic Bulge 

The thickness of the Galactic bulge (ratio of minor to major axis) is a very sensitive 
probe of the environmental conditions in which the Milky Way's bar formed and evolved, 
and disagreements as to the thickness have been a catalyst to significant disagreements as 



to the nature of the Milky Way's bar (Lopez-Corredoira et al. 2011). Further, detailed 



investigations of the Galactic bulge luminosity function require not only the first moment of 
the distance distribution, but also the second moment as well. 



a 



We estimate the geometrical thickness of the Galactic bulge in units of distance modulus, 
as follows: 



cr, 



a 



l,RC 



I,RCfl 



a 



{V-I),RCfl I ' 



where (Ti^rc is the measured brightness dispersion of the RC, ai^Rc,i 
intrinsic magnitude dispersion of the RC in /, and -Rf x (a 



{v-i) ^(y-i),RCfl 



extinction component, with Ri = Aj/E{V — I) measured as per Equation 28 



(33) 

0.09 is the estimated 
is the differential 
Due to 



the fact the measurement error on aj^jic is large relative to the variations, we smooth the 
measurements by replacing each value of cr^ with the mean of all the values of cr^ located 
within 20'. The values of function of direction are plotted in the bottom panel of 

Figure 20 We also show the scatter of cr^ for latitudinal and longitudinal stripes on the 



bottom-left and bottom-right panels of Figure 21, respectively. The bottom-left panel shows 



that for sightlines close to the plane, is minimized near / = 0°, the expected behavior 



for a triaxial ellipsoid 

near / = ±5°. Nataf et al. (2010) used the increase in aj^RC toward (/ 



Conversely, for sightlines further from the plane, is minimized 

0°, \b\ > 4.75°) to 



infer the existence of the double RC. This increased geometric dispersion with separation 
from the plane for sightlines along the minor axis is matched by N-body models with an 
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X-shaped bulge ( |Li k Shen|[2012[ |Ness et al.|r2012D . Further from the minor axis, we suspect 
that the high values of for / > 8° are due to disk contamination. The bottom-right panel 
of Figure 21 shows that o"^ is correlated with separation from the plane along the minor 
axis (as would be expected of a bulge), and becomes anti- correlated with separation from 
the plane at large separations from the minor axis. The latter is the expected behavior if 
disk contamination is high for / = 10°: sightlines closer to the plane will probe more of the 
disk and thus a larger range of distances. We note that the very low values of close to 
(Z, h) = (0°, —2°) are more sensitive to possible errors in our zero-point calibrations. 



Table 2: Values of the mean distance modulus (assuming Mj^ptc = —0.12), distance modulus 
dispersion, reddening, differential reddening, and reddening law (Equation 28) for four 



sightlines that are the subject of intensive HST observations ( Sahu et al. 2006 Brown et al. 
2010 Clarkson et al. 2011). E(y — /) and aE{v-i) for these sightlines were measured using 



circles centered on these field centers with radii of 3'. 



Field Name 


I (deg) 


b (deg) 






E{V - I) 


(^E{V-I) 


RjKVI 


Stanek's Window 


+0.25 


-2.15 


14.53 


0.17 


1.04 


0.08 


0.375 


SWEEPS 


+ 1.26 


-2.65 


14.52 


0.20 


0.79 


0.10 


0.362 


Baade's Window 


+ 1.06 


-3.81 


14.54 


0.24 


0.67 


0.04 


0.351 


OGLE 29 


-6.75 


-4.72 


14.76 


0.29 


0.67 


0.00 


0.321 



The spatial distribution of shows evidence for at least three structures. The inner 
part of the Galactic bulge probed by the OGLE-III survey, for which (|/| < 4°, |6| < 3°), is 
shaded magenta in both panels of Figure 20, and appears consistent with a perfect triaxial 
ellipsoid. Further from the plane, but within \l\ < 5°, we see signs of an X-shaped structure, 
with the geometric dispersion being maximized along the minor axis. For sightlines more 
than 5 degrees away from the minor axis, the geometric dispersion is measured to increase 
rapidly, which may be due to a much higher disk contamination fraction. The geometric 
dispersion toward four bulge fields of high scientific interest, including Baade's window and 
Stanek's window (Stanek 1998), are listed along with other parameters in Table [2] Finally, 
we note that the geometric dispersions measured here apply only to the stellar population 



adequately traced by the RC. Metal-poor bulge stars are kinematically hotter ( 


Babusiaux et 


al. 


2010 


De Propris et al. 


2011 


Pietrukowicz et al. 


2012 


Uttenthaler et al. 


2012 


), and will 



thus have a larger geometric dispersion. 
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ls. Summary of Data 

The results of this work are available for download on the OGLE webpage We briefly 
summarize the format here. The flrst table is intended for observers, of which we show 
a cropped version as Table |3| It includes values of the reddening, extinction, differential 
reddening, reddening law, mean distance modulus, and distance modulus dispersion for each 
sightline. The OGLE webpage includes a GUI to facilitate retrieval of these values for a 
given coordinate. The determinations of /i, o"^, Nrc and full error matrix thereof are also 
available in a separate table. We show a cropped version as Table |4j 

Table 3: Coordinates, extinction, reddening, differential reddening, mean distance mod- 
ulus, distance modulus dispersion, reddening law and Quahty Flag for each of ~9,000 
sightlines studied in this work. A 20' smoothing is applied to the values of the mean and 
dispersion of the distance modulus distribution. A Flag of "0" means a reliable measurement. 
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E{V - I) 


(yE(v~i) 


^^ 
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Flag 


-10.06 


-4.32 


0.89 


0.78 


0.07 




17 


0.34 


0.29 





-9.94 


-4.24 


1.02 


0.87 


0.09 




19 


0.39 


0.32 





-9.91 


-4.57 


0.85 


0.72 





lAl 


16 


0.35 


0.31 






Table 4: Structural parameters for sightlines deemed reliable, with full error matrix. 
The symbol "S"' denotes the l-cr measurement error on the variable, as opposed to the 
more standard "cr" to avoid confusion with a^, and C denotes the correlation. Due to 
measurement errors, approximately 2% of sightlines have a best-flt negative variance to the 
distance modulus distribution. 
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0.13 
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0.07 


0.04 
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0.14 
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14.84 


0.13 
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0.07 


0.05 


54 


0.22 


0.23 


0.72 
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14.92 


0.20 


394 


0.07 


0.06 


74 


0.33 


0.37 


0.82 



ht tp : / / ogle . astrouw . edu .pi/ 



- 51 - 



14. Discussion and Conclusion 



Our solution to the observational problem of the non-standard VI extinction toward 
the inner Galaxy mitigates what has been one of the dominant sources of uncertainty in 
studies of the Galactic bulge. The extinction law is on average steeper in the optical with 
~30% variations superimposed, and is well- fit by the relation Aj = 0.74:65xE{V — I) + 
1.3700xE{J — Ks). The residuals to the extinction fit is now reduced to no more than 0.06 
mag, and the estimate of Rq = 8.20 kpc is consistent with there being no bias in our fit to 
the extinction law. In the course of making these measurements, we have also measured that 
differential reddening averages (~9% of total reddening for small fields), and the intrinsic 
luminosity parameters for the bulge RC. These will be of use to future bulge studies. 

The mean value of Ai/E(y — I) = 1.217 suggests Ry = 2.5, the mean value of Rjkvi = 
0.3433 suggests Ry = 2.6, the mean value of Rvi,BVsid = 1-179 suggests Ry = 2.5, and 
thus all the measurements investigated here are consistent with an Ry ^ 2.5 extinction 
curve. Measuring the extinction curve in other bandpasses could potentially have major 
implications for cosmology. Our inferred extinction curve is consistent with the values of 
Ry ~ 2.5 inferred in studies of the extinction toward extragalactic type la SNe by |Guy et 



al. (2010) and Chotard et al. (2011). In particular, the hierarchical Bayesian analysis of 
Mandel et al. (2011) found that the extinction law toward SNe la went as Ry = 2.5 — 2.9 



for Ay < 0.4, and steepened at higher extinctions. Falco et al. (1999) also reported a 



range in the extinction laws of 23 lensed galaxies of 1.5 < Ry < 7.2. Since none of the 



combined 51 measurements of Cardelli et al. (1989) and O'Donnell (1994) reach such low 
values of Ry, there is no reason to expect that extrapolation of this empirical law will 
behave adequately in domains that lie well beyond its calibration. Indeed, though the values 
of Ai/E{V — I) = 1.217 and Rjkvi = 0.3433 are consistent with each other, the change in 



Ai/E{V — I) as Rjkvi changes does not go at the rate predicted by Cardelli et al. (1989) 



and O'Donnell (1994), as shown in Figure 13 This domain of the extinction law therefore 



warrants further investigation. Some constraints could be extracted by combining our results 



with the recent study of Nidever et al. (2012), who measured reddening values for the color 
([3.6/x] - [4.5/i]). 

Our measurements of the number counts, brightness dispersion, mean brightness and 
full error matrix thereof for ~9,000 RC centroids toward the bulge may be one of the most 
potent means for constraining the structural parameters of the Galactic bulge. In Sections 



10 and 11 , we have sketched how these data could be used to constrain the morphology and 



mass of the bar, without going to the full formalism employed by Dwek et al. (1995), Stanek 



et al. (1997) and Rattenbury et al. (2007a). Moreover, as the data have improved, it is now 



time for the models to improve as well. The use of N-body models by Rattenbury et al. 



( |2007bD , |Shen et al.| ( [20lo| , [Martinez- Valpuesta fc Gerhard] ( [20ll| , and [Ness et al.| ( [20l2| ) 

are encouraging steps in that direction. 
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